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^ : Abstract 

This paper presents an analysis of a model problem, consisting of two inter- 

acting rigid rings, for the rotation of molecules in liquid '^He. Due to Bose 

O ] symmetry, the excitation of the rotor corresponding to a ring of N helium 

[ atoms is restricted to states with integer multiples of N quanta of angular 

'c/2 ' 

momentum. This minimal model shares many of the same features of the 
rotational spectra that have been observed for molecules in nanodroplets of 
~ 10'^ — 10^ helium atoms. In particular, this model predicts, for the first 
' time, the very large enhancement of the centrifugal distortion constants that 

' have been observed experimentally. It also illustrates the different effects of 

' increasing rotational velocity by increases in angular momentum quantum 

Q ' number or by increasing the rotational constant of the molecular rotor. It 

is found that fixed node, diffusion Monte Carlo and a hydrodynamic model 
provide upper and lower bounds on the size of the effective rotational constant 



X 
J3 



of the molecular rotor when coupled to the helium. 
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The spectroscopy of atoms and molecules dissolved in helium nanodroplets is a topic of 
intense current interest One particular, almost unique feature of this spectroscopic 

host is that even heavy and very anisotropic molecules and complexes give spectra with 
rotationally resolved structure [Q. This spectral structure typically corresponds to thermal 
equilibrium, with T ^ 0.38 K, and has the same symmetry as that of the same species in the 
gas phase The rotational constants, however, are generally reduced by a factor of up 

to four or five, while the centrifugal distortion constants are four orders of magnitude larger 
than for the gas phase . These large changes clearly reflect dynamical coupling between 
the molecular rotation and helium motion. At present, there are at least four different models 
proposed for the increased effective moments of inertia, at least two of which have reported 



quantitative agreement with experiment [H-ll]. The large observed distortion constants have 



not yet been quantitatively explained, and the most careful attempt to date to calculate them 
(for OCS in helium) gave an estimate ~ 30 times smaller than the experimental value 

The highly quantum many body dynamics of this condensed phase system has made it 
difficult to achieve a qualitative understanding of the observed effects. In cases like these, 
simple models can provide insight, especially if the lessons learned can be tested against 
more computationally demanding simulations that seek, however, to provide a first principles 
treatment of the properties of the system of interest. In this paper, one very simple model 
system will be explored that seeks to model the coupling of a molecular rotor to a first 
solvation shell of helium. The existing models for the reduced rotational constants agree 
that most of the observed effect comes from motion of helium in the first solvation shell. 
Some of the qualitative features of this model were discussed previously , but quanitative 
details were not persued in that work. 

The 'toy' model considered consists of a planar rotor coupled to a symmetric planar 
ring of helium atoms. This model problem can be solved exactly, and can reproduce the 
size of the observed reductions in the rotational constant AND the size of the centrifugal 
distortion constants. This is the first time, to the authors knowledge, that the large effective 
distortion constants of molecules in liquid helium has been reproduced. Further, this model 



clearly resolves a confusion about the sign of the centrifugal distortion constant. Based 
upon the expected decreased following of the helium with increasing rotational angular 



rigid rotor, ie., that the effective centrifugal distortion constant should be negative, in conflict 
with experimental observations. The present model demonstrates, however, that opposite 
behavior is expected when the rotational velocity of the rotor is increased by increasing the 
rotational quantum number (where an increased angular anisotropy and following of the 
helium is predicted) or when the rotational constant of the isolated rotor is increased (where 
decreased angular anisotropy and following of the helium is predicted). The present model, 
therefore, rationalizes both the observed depenence of the increased moments of inertia 
on the rotational constant of the isolated molecule and the observed centrifugal distortion 
constants. 



We will consider a highly abstracted model for rotation of a molecule in liquid helium. 
The molecule will be treated as a rigid, planar rotor with moment of inertia Ji. The ori- 
entation of the molecule is given by 6i The liquid helium is treated as a ring of helium 
atoms that forms another rigid, planar rotor with moment of inertia I2 and with orientation 
given by 62- Because of the Bose symmetry of the helium, he helium rotor can only be 
excited to states with A^^ units of angular momentum. The lowest order symmetry allowed 
coupling between the molecule and the helium ring is given by a potential V cos [N{6i — 62)]. 
Any coupling spectral components that are not multiples of A^ will lead to mixing of states 
that are not allowed by Bose symmetry, which is forbidden in quantum mechanics. The 
Hamiltonian is given by: 



We define B12 = the rotational constants for the uncoupled rotors. We can separate 
the above H by introducing the two new coordinates: 



velocity one can argue that the rotational spacing should increase faster than for a 



I. THE TOY MODEL 




(1) 



3 



- _ Ii9i + I2O2 
h + h 



(2) 



in which we have: 



2(/i + h) de\ 



2 + 




l_\ ^ 



+ Vcos{N9) 



(3) 



9 is the variable conjugate to the total angular momentum; ^ is a vibrational coordinate. 
We define -Brigid — 2(h+i2) -^""^i ~ ¥ !^)" eigenstates of H separate into a 
product: 



7j;{9,9)^e''%{9) 
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J is the quantum number for total angular momentum. It would appear from the separable 
H that the energy could be written as an uncoupled sum of a rigid rotor energy, -Brigid ( 
not J(J+ 1) because we have a planar rotor), and a 'vibrational' energy that is independent 
of J. However, the energies are not simply additive, due to the fact that the boundary 
condition for 9 is J-dependent. When ^2 alone is changed by any multiple of 2t^/N, t/j must 
be unchanged. However, a change of —2tt/N in 92 results in a change of —{2t:/N) {I2/ {h+h) 
in 9 and +27r/N in 9. Thus, the Bose symmetry of the helium ring is satisfied by taking as 
the boundary condition for ipy-. 



As a result, the 'vibrational energies' and eigenfunctions are a function of the total angular 
momentum quantum number, J. Note that the boundary condition is periodic in J, with 
period given by A^(/i + h)/h- 

The J dependence of the boundary condition of the vibrational function is rather unfa- 
miliar in molecular physics. This dependence can be removed by a Unitary transformation 
of the wavef unction: 




(5) 
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The boundary condition on the transformed function is iplmv (f + if) — i^'rmv (^) • The 
transformed hamiltonian, H\ is given by: 

H' = ^^{J-Lf + ^L' + Vcos {NO) (7) 

where J = ih-^ is the operator for total angular momentum, and L = is the angular 
momentum of the helium ring relative to a frame moving with the rotor. This form of 
the Hamiltonian closely resembles those widely used in treatment of weakly bound com- 
plexes Expansion of the (J — L)^ gives a Coriolis term that couples the overall rotation 



to the vibrational motion, which makes the nonseperability of these motions evident. This 
form, however, hides the periodicity in J of this the coupling. 

We will now consider the two limiting cases. The first to consider is that where |V^| S> 
A^^^^^j. In this case, the potential can be considered harmonic around 6 = 2TTk/N,k = 
. . . N — 1 when V is negative (and shifted hj tt/N for V positive) with harmonic frequency 
v = 2^y''2A^^|l^| (^j- + -^). The wavefunction will decay to nearly zero at the maxima of 
the potential, and changes in the phase of the periodic boundary condition at this point 
(which happens with changes in J) will not significantly affect the energy. In this limit, the 
total energy is E{J,v) = -Brigid<^^ + hui^v + 1/2), and we have a rigid rotor spectrum with 
effective moment of inertia Ji + /2. While there are N equivalent minima, Bose symmetry 
assures that only one linear combination of the states localized in each well (the totally 
symmetric combination for J = 0) is allowed, and thus there are no small tunneling splittings, 
even in the high barrier limit. The total angular momentum is partitioned between the 
two rotors in proportion to their moments of inertia, i.e. < Ji >= hJ ■ /i/(/i + I2) and 

< J2 >=hj -hlih + h). 

We will now consider the opposite, or uncoupled rotor, limit. The eigenenergies in this 
case are trivially £{1711,1712) = Bim\ + B2m\ with eigenf unctions ip = exp {irriiOi + 1771262)- 
mi can be any integer, while ni2 = Nk, where k is any integer. Introducing the total angular 
momentum quantum number J = mi + m2, we have E{J, 7712) = Bi{J — 7/12)^ + B2 m\. The 
lowest state for each J has quantum numbers mi = J — Nk and m2 = Nk, where k is the 



nearest integer to BiJ /N{Bi + B2). Treating the quantum numbers as continuous, we have 
E = i^rigidJ^, e.g., the same as for rigid rotation of the rotors. However when we restrict J 
to integer values, for J < N{Bi + B2) / {2Bi), the energy spacing will be exactly that of a 
rigid rotor with rotational constant Bi. In general, as a function of J, the uncoupled ground 
state solutions follow the rigid rotor spectrum in but with a series of equally spaced 
curve crossings when the lowest energy 1112 value increases by as J is increased by one 
quantum. These crossings allow the total energy to oscillate around that predicted for a 
rigid rotor with moment of inertia Ji + l2- 

II. NUMERICAL RESULTS 

Having handled the limiting cases, we can now turn our attention to the far more in- 
teresting question, which is how does the energy and eigenstate properties change as V is 
continuously varied between these limits. We note that changing the sign of V is equivalent 
to translation of the solution by A^^ = tt, and thus we will only consider positive values 
of V explicitly. We also note that the eigenstates are not changed, and the eigenenergies 
scale linearly if Bi,B2, and V are multiplied by a constant factor. As a result, we will 
take i?2 = 1 to normalize the energy scale. The solutions for finite values of V were cal- 
culated using the uncoupled basis and the form of H given in Eq. |l| with fixed values of 
mi + m2 = J and m2 = Nk. For each value of J, the matrix representation for if is a 
tridiagonal matrix, with diagonal elements given by the energies for the uncoupled limit, 
and with off-diagonal elements given by 1^/2. Numerical calculations were done using a finite 
basis with k = —15, —14, ... 15. 

Using Bi = B2 and = 8, we have calculated the lowest eigenvalues of H for J = 0, 1, 2 
and used these, by fitting to the expression E{v, J) = Eoiy) + BcsJ^ — DcsJ"^, to determine 
Bcs and D^s- Figure ^ shows the value of i?efr as a function of V (both in units of B2). It 
can be seen that Bcs varies smoothly from Bi to -Brigid with increasing V, and is reaches a 
value half way between these limits for V ~ A^i?2. 
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In order to rationalize this observation, we will now consider a Quantum Hydrodynamic 
treatment for the rotation |T5|. Let the ground state density be p{6) = \iIj^{6)\'^. Now let the 



molecule classically rotate with angular velocity uj. To first order in uj, p will not change (i.e. 
we will have adiabatic following of the helium density for classical infinitesimal rotation of 
the molecule). However, the vibrational wavefunction, ip^ will no longer be real, but instead 
will have an angle-dependent phase factor whose gradient will give a hydrodynamic velocity. 
Solving the equation of continuity: 

where r is the radius of the helium ring, gives solutions of the form: 

(l__£_) (9) 

where C is an integration constant. We determine C by minimizing the kinetic energy 
averaged over 9. This gives: 

27r 

C = (10) 

and a kinetic energy: 

In the case of a uniform density, p = (27r)^^ and AEk = 0. As the density gets more 
anisotropic, the integral becomes larger and AEk becomes larger, approaching the value 
for rigid rotation of the helium ring when p has a node in its angular range. We define 
the hydrodynamic contribution to the increase in the moment of inertia of the heavy rotor 
due to partial rotation of the light rotor by AEk = ^AIy^uj"^. It is interesting to note that 
for the above lowest energy value of C, we have Jq^ v{6)d6 = 0. {i.e. that the solution is 
'ir rotational') and that the net angular momentum induced in the helium is A/ho;. The 
lowest energy solution of the three dimensional Quantum hydrodynamic model satifies these 



conditions as well ||TO| , P^ 



The hydrodynamic model can be tested against the exact quantum solutions. Define 
A/eff as the effective moment of inertia for rotation (as calculated from -Bes) minus the 
moment of inertia for the molecular rotor. A/off will grow from for uncoupled rotors 
to I2 as the coupling approachs the rigid coupling limit of high V . In the hydrodynamic 
model, A/gff = AJh. Figure g shows a plot that compares A Jeff and AJh as a function of V . 
Each has been normalized by l2- They are found to be in qualitative agreement for the full 
range of V , though the exact quantum solution is systematically below the hydrodynamic 
prediction. We note, however, that for the assumed parameters, the speed of the molecular 
rotor is equal to that of the helium rotor, while the hydrodynamic treatment assumed a 
classical, infinitesimal rotation of the molecular rotor. 

The size of A/og is determined by the degree of anisotropy of the ground state density in 
the vibrational displacement coordinate 9. If Ji is decreased at fixed I2 and V , the effective 
mass for ^, which is (Jf^ + l2^)^^ will also decrease, which will decrease the anisotropy 
produced by V . Fig |^ shows how the normalized A/gg and A/h vary as the molecular 
rotational constant, -Bi, changes from to 2B2. This calculation was done for V = 100, 
close to the value corresponding to maximum difference of AI^s and AJh for Bi = B2. 
This plot demonstrates that the hydrodynamic prediction becomes exact in the limit that 
Bi 0, i.e., in the case that the assumption of infinitesimal rotational velocity of the 
molecule holds. However, it substantially overestimates the increase effective moment of 
inertia when Bi > B2. This decrease in the increase moment of inertia with increasing 
rotational constant of the heavy rotor is the effect previously interpreted as the breakdown 
of adiabatic following in the literature on the rotational spectrum of molecules in liquid 



Helium 12,10,13 



Figure ^ shows a plot of D^s as a function of V for Bi = B2 = 1. D^s = is zero in both 
limits, and has a maximum value near the value of V at which B^s is changing most rapidly. 
It is interesting to explicitly point out that this Dcs value arises entirely from changes in the 
angular anisotropy of the helium density with J, as the model does not allow for an increase 
in the radial distance of the helium, which has previously been considered 0. Further, the 
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peak value of Dea ~ 1.8 • 10~'^-Bi is in remarkably good agreement with the ratio of 

to the gas phase molecular rotational constant observed for a number of molecules in liquid 

helium. For example, for OCS this ratio is found to be 2 • 10^"^ while for HCCCN, the 



same ratio was found to be 1 ■ 10 ^ 



We can gain further insight by examining the rotational energy systematically as a func- 
tion of J. Figure |^ shows the rotational excitation energy (-E'(0, J) — E{0,0)) divided by 

as a function of J. The calculations were done with V = 100. The rotational excitation 
energy approaches that of the -Brigid</^ for high J. Further, it reaches this value for J equal 
to multiples of A^(/i + I2)/ 121 which matches the periodicity of the boundary conditions for 
ip^. J values that lead to the same boundary conditions for ip^ will differ in energy only by 
the eigenvalues of i^r, and thus it follows from Eq. ^ that of a rigid rotor with rotational 
constant -Brigid- For the first half of each period in J, ip^ is found to increase in its anisotropy, 
and therefore the energy increases, as J is increased (See Fig. |]). This can be understood 
when one considers the fact that for J = N{Ii + I2)/ {212)1 the boundary condition is that 
iI)^{2t{/N) = — '(/'v(O), i.e. the wavefunction will be real but have nodes in the interval 
[0,27r]. 

Classically, the molecular rotor is characterized by its rotational angular velocity, u = 
2BiJ. However, we see that the quantum treatment of the two coupled rotors gives opposite 
results when u is increased by increasing either Bi or J. For increases in Bi, the 'degree 
of following' of the light rotor decreases for fixed potential coupling, as seems intuitively 
reasonable. However, for increases in J, the anisotropy of the potential and thus the 'degree 
of following' initially increases, and thus so does the effective moment of inertia of the coupled 
system. This behavior continues until one passes through a resonance condition where the 
helium can be excited by transfer of quantum of angular momentum from molecular 
rotor to the helium. This resonance condition is missing from the classical treatment of the 
coupling between the rotors, where the angular velocity of the molecular rotor is treated as 
a fixed quantity, u, which is one of the parameters of the problem. 
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III. NODAL PROPERTIES OF SOLUTIONS 



It is possible to calculate the rotational excitation energies of clusters of helium around 
a molecule by use of the Fixed Frame Diffusion Monte Carlo (FFDMC) method |jl2|. As 



in most DMC methods, this method should yield (except for statistical fluctuations) a 
upper bound on the true energy, finding the optimal wavefunction consistent with the nodal 
properties that are imposed on the wavefunction by construction. In the case of FFDMC, 
the nodal planes are determined by the free rotor rotational wavefunction for the molecule 
alone, i.e. that the sign of the wavefunction (which is taken to be real) for any point in 
configuration space is the same as that of the rotor wavefunction at the same Euler angles. 

We can examine the exact solutions of our toy problem to gain insight into the accuracy 
of the nodal planes assumed in FFDMC. The wavefunctions we have considered up to now 
are complex, but because of time reversal symmetry, the solutions with J and — J rotational 
quantum numbers must be degenerate. Symmetric combination of these solutions just gives 
the Real part of J solution, and the antisymmetric combination the Imaginary part. The 
real part is given by: 

^/^^(^i, 02) = cos{J9i)Re{ij'{92 - 9i)) - sin( J^2)/m(^'(^2 - ^i)) (12) 

where 

ilj'iAO) =Y.4 eM^kNAO) (13) 

k 

and c{ are the eigenvector coefficients obtained from diagonalization of the real Hamiltonian 
matrix in the uncoupled basis. Examination of the numerical solutions reveals that for 
J < N, Re{ip') has no nodes, while Im{ip') has nodes at equal to integer multiples of 
tt/N. Thus, if Im{ilj') = 0, then the solution would satisfy the FFDMC nodal properties 
exactly. However, for finite Im{ilj'), the nodal surfaces, rather than being on the planes 
6i = constant, are modulated N times per cycle along the di = constant line. For = 80 
and Bi = B2 = 1, the maximum value of Im{ip') is about 4% of Re{ip'), and growing 
approximately linearly for low J. 
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In order to test the quantiative implications of this error in the nodal properties, impor- 
tance sampled DMC calculations have been done for the present two rotor problem. The 
explicit DMC algorithm given by Reynolds et al. |T^ was used with minor change ||19| . The 



guiding function, -j/^t, which determines the nodes, was selected as cos( J6'i)?/'°(6'i — ^^2), where 
is the real, positive definite eigenstate for the J = problem. The rotational constant, 
i?DMC is defined as the DMC estimated energy for J = 1 less the exact ground state eigenen- 
ergy for J = 0, and will be (except for sampling and finite time step bias) an upper bound 
on the true B value calculated earlier. The points plotted in figure |l| are the calculated 
values of -Bdmc with the estimated 2a error estimates. It is seen that the fixed node DMC 
estimates of B^s are excellent for low values of V , but underestimate the contribution of the 
Helium ring to the effective moment of inertia as it is coupled more strongly to the rotor. 

IV. RELATIONSHIP WITH A MORE REALISTIC MODEL 

In a series of insightful lectures, Anthony Leggett analyzed the properties of the ground 



state of N Helium atoms confined to an annulus of radius R and spacing <^ -R p5 
The walls of the annulus are allowed to classically rotate with angular velocity uj. While 
not stated explicitly, the walls of the annulus couple to the helium via a time dependent 
potential, which is static in the rotating frame. As such, our rotating diatomic molecule 
can be considered as a special case of the problem treated by Leggett. If one transforms to 
the rotating frame, the quantum hamiltonian is the same as for the static [u = 0) problem. 
However, the boundary condition for the wavefunction in this frame is given by [|l^[Eq. 
(2.10)]: 

^o(0i, ^2 . . . % + 27r . . . 0jv) = exp {2mmR^uj/h) ^'^{61, 62... Oj... On) (14) 

In making a comparison to the results of the toy model, we note that for this system I2 = 
NmR'^ (the classical moment of inertia for the helium) and J = u!{Ii + 12)/^. Substitution 
shows that the phase factor in Eq. is identical to that derived above for Eq. ^ Note, 
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however, that Eq. refers to moving one hehum atom by 27r, while Eq. ^ refers to motion 
of all N helium atoms by 27t/N. Motion of all helium atoms by 2tt will result in a phase 
factor of 2TiiN mR'^ uj / h = 2Ttil2J / {Ii + I2) in both treatments. 

Leggett considered the change in helium energy produced by rotation of the walls. Let 
Eq be the ground state energy for the static problem, and Eq{uj) the ground state energy 
in the rotating frame. The ground state energy in the laboratory frame is given by [Eq. 
(2.12)]: 

^lab = ^0 + ^/2^'-[^0M-^0] (15) 

For the ground state of Bosons, we further have that Eq{uj) > Eq, with equality only when 
u equals integer multiples of uq = h/mR'^ since the nodeless state has the lowest possible 
energy. At u = kujQ, the helium rigidly rotates with the walls. This agrees exactly with 
the numerical results of the toy model, as shown in Figure In making comparisons with 
this model, one should remember that Eig^y, does not include the kinetic energy of the walls 
(rotor). Thus the more general treatment of Leggett supports one of the central insights of 
the toy model, that the large effective distortion constants for molecular rotors in helium is 
a consequence of an increased helium following of the rotor with increasing angular velocity, 
which in turn is a direct consequence of the u dependence of the single- valuedness boundary 
condition in the rotating frame. 

The moment of inertia for the ground state of the helium can be defined by: 

Leggett defined the "normal fraction" of the helium by the ratio I/I2, which is equal to unity 
if Eq{(jj) is independent of as ^ 0. This will occur if the wavefunction has 'nontrivial' 
nodal planes, since the phase of the wavefunction can be changed discontinuously at a 
node without cost of energy. Nodal plans associated with overlap of particles, however, are 
'trivial' in that the phase relationship on each side of the node is determined by the exchange 
symmetry of the wavefunction, and thus cannot be used to match the boundary conditions 
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without extra cost of energy. In our toy problem, when V is very large, the vibrational 
wavefunction becomes localized, introducing near nodes at the maxima of the potential, and 
as a result the ground state is described by a near unity normal fraction; we have what 
Leggett refers to as a 'normal solid'. Conversely, as the uncoupled limit is approached, the 
helium ring does not contribute to the kinetic energy of the lowest rotational states and 
we have / ^ 0, and we have zero normal fraction (i.e. the helium has unity superfluid 
fraction). Following Leggett's definition, one finds that the normal fraction is given by 
A/eff//2. Thus, Figure ^ can thus be interpreted as the normal fluid fraction for the ground 
state as a function of the strength of the potential coupling. Leggett's analysis is based upon 
a classical treatment of the motion of the walls, which implies Ii ^ I2, in which limit the 
hydrodynamic model exactly predicts the normal fluid fraction. 
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FIGURES 



Effective Rotational constant for two rotors. B =B =1., N=8. 
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FIG. 1. Effective Rotational Constant for two coupled rotors as a function of the interaction 
potential strength. Bi = B2 = 1, and the second rotor can only be excited to states with multiples 
of 8 quanta. The individual points are the B effective values calclated by a fixed node, Diffusion 
Monte Carlo calculation. The error bars on these points are the estimated 2a sampling error. 
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FIG. 2. Increase in effective moment of inertia of molecule, A^s due to coupling to rotor made 
of 8 helium atoms. A/hydro is the same quantity estimated by the hydrodynamic model. Both are 
calculated as a function of the potential coupling strength, and the results normalized to the rigid 
rotor moment of inertia of the 8 helium rotor 
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Change in moments of inertia with molecular rotational constant 
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FIG. 3. Same as Figure |2|, except as a function of the rotational constant of the molecule, 
normalized to the rotational constant of the 8 helium rotor. 
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Effective Distortion Constant for two rotors. B = B = 1 , N = 8 
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FIG. 4. The effective centrifugal distortion constant, -Dcflf) for molecule coupled to ring of 8 
helium atoms as a function of the strength of the coupling, V. Both Dcff and V are normalized to 
the rotational constants of the molecule and 8 helium rotor, which are taken as equal. 
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FIG. 5. The rotational excitation energy, AE, divided by as a function of the total rotational 
angular momentum quantum number, J. Calculated with Bi = B2 = 1 and V = 100. With rigid 
following of the helium, the plotted quantity should equal -Bi-igid? which is indicated in the figure. 
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Changes in wavefunction with rotational excitation 
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FIG. 6. The absolute value of the vibrational wavefunction as a function of the relative ori- 
entation between molecule and 8 helium rotor. This quantity is plotted for angular momentum 
J = 0,1... 8. 
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